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Abstract 

Filter banks and wavelet decompositions that employ recursive filters have been considered 
previously and are recognized for their efficiency in partitioning the frequency spectrum. This 
paper presents an analysis of a new infinite impulse response (HR) filter bank in which these 
computationally efficient filters may be changed adaptively in response to the input. The fil- 
ter bank is presented and discussed in the context of finite-support signals with the intended 
application in subband image coding. 

In the absence of quantization errors, exact reconstruction can be achieved and by the 
proper choice of an adaptation scheme, it is shown that IIR time-varying filter banks can yield 
improvement over conventional ones. 


1 Introduction 

Subband image coding is a well-known technique in which an input is split into a small number of 
subband images, each of which is associated with a different region in the spatial frequency plane. 

The subband images are decimated to their Nyquist rates prior to being quantized and coded. 

In reconstruction the images are decoded, upsampled, and then merged together. The process 
in which the input is split and merged is called analysis/synthesis. There are many performance 
and systems related design issues associated with analysis/synthesis as discussed in [2], [4], [7] 
and elsewhere. Recursive analysis and synthesis filter banks have been shown, to be particularly 
attractive because they can be designed to reconstruct the input exactly in the absence of coding 
and can achieve tremendous computational efficiency relative to comparable FIR systems [4], 

[5]. These filter banks are typically two-band systems that serve as building blocks for uniform 
and non-uniform band tree structured systems. However, parallel form recursive filter banks 
have also been shown to work well [1 1] i In such a case, the subband images are complex valued. 

1 This work was supported in part by the National Science Foundation under contract MIP-9116113 and the 
National Aeronautics and Space Administration. 
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Thus to avoid increasing the data rate, Hermitian symmetry is exploited to counter the sample 
increase due to the presence of complex numbers in the subband images. 

When subbands are quantized, distortion results and the degree of the distortion is inversely 
proportioned to the bit rate. The quantization noise generates three types of interrelated dis- 
tortions: 

1. Aliasing. This is due to the partial loss of aliasing cancellation in the synthesis section that 
arises from quantization. It appears subjectively as blurring in the reconstructed images. 

2. Ringing distortion. This appears as amplitude oscillations (or ringing) in the vicinity of 
edges in the image. As discussed in [4], its cause may be directly related to the ripples in 
the step response characteristics of the analysis/synthesis filters. 

3. In-band spectral distortion. This distortion is attributed to the deviation in the spectral 
magnitude and phase of subband images excluding the inter-band effects of aliasing and 
ringing. In essence, this form of distortion accounts for all remaining deviations when 
aliasing and ringing effects are removed. 

Much attention has been given to minimizing the perceptually important coding distortions 
in image coding systems. For low bit rate subband image coding systems, aliasing and ringing 
distortions are often cited as being objectionable. Efforts to reduce these effects have been 
considered previously but with limited success [4], [5], [12]. The difficulty is that the ripples 
of the step response of the filters, which are the source of the ringing distortion, are needed to 
achieve good filter magnitude response characteristics. Good spectral magnitude characteristics 
are important for reducing the effects of aliasing in the presence of coding. Aliasing distortion 
appears as blurring in the image and is subjectively noticeable. If filters with good magnitude 
characteristics are used to reduce the aliasing, then ringing distortion results that is noticeable 
and objectionable. As to which is worse is a matter of opinion. It is clear, however, that the 
elimination or reduction of both is desirable. This dilemma is a consequence of the Gibbs Phe- 
nomenon which implies that is it impossible to have filter banks with good (i.e. monotonic) step 
response characteristics and good magnitude response characteristics within the conventional 
paradigm. However, it is possible to overcome this seemingly fundamental difficulty by using 
time varying filter banks as we shall see. The first introduction of the concept of time-varying 
filter banks appears in [3] and is based on the time domain theory and design methodology 
presented in [15], [14], and [13]. The approach taken in [3] has some attractive features: time- 
varying capability is conceptually simple in this framework and system delay can be controlled 
very easily. On the other side of the coin, the approach is inherently limited to FIR systems 
and the system design task can be difficult for practical filters. 

In this paper, we introduce the concept of time- varying recursive filter banks and demonstrate 
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Figure 1: Two-Band Analysis/Synthesis System 

that they provide the additional degree of freedom that allows them to overcome the step- 
response-magnitude-response constraint associated with the Gibbs Phenomenon. We derive the 
conditions under which we can switch between filters of desirable properties. Furthermore, 
this new time-varying recursive filter preserves the exact reconstruction property and maintains 
very high computational efficiency. We presented this idea at the 1992 DSP Workshop [16] and 
discussed some of the analysis issues at the 1992 Asilomar Conference [17]; both [16] and [17] 
are rather limited due to page length constraints. 2 

This paper begins with a discussion of the exact reconstruction IIR filter bank after which 
the time-varying conditions are derived. The equations for time-varying analysis/synthesis are 
derived and a design strategy is developed that allows for much improved performance over 
conventional filter banks. 

2 Exact Reconstruction IIR Filter Banks 

The starting point for this work is the two-band analysis/synthesis system shown in Figure 1 
where {//oC*), H\ (z)} and {G'o(z), G\{z)} are the pairs of lowpass/highpass analysis/synthesis 
filters. The analysis filters have the form: 


Ho(z) 

= PoC^ + z-'fMz 2 ) 

(1) 

Hi(z) 

= P 0 (z 2 ) - z- 1 p^z 2 ) 

(2) 


and may be implemented using the well-known polyphase structure shown in Figure 2 where 
P 0 (z) and P\{z) are the analysis polyphase filters and Qa(z) and Qi(z) are the synthesis 
polyphase filters. 

The exact reconstruction property of this recursive filter bank is evident from an examination 
of the polyphase structure. Observe that the criss-cross networks in the analysis and synthesis 
sections of Figure 2 are simply two-point DFT butterflies and that the cascade of the two 


2 The authors recently became aware of similar work developed independently in Norway by Husoy and Aase [18]. 
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Figure 3: Two-Band Polyphase Implementation with DFT Butterflies Omitted 

results in an identity system. Therefore the combined analysis/synthesis polyphase system can 
be simplified as shown in Figure 3. By inspecting Figure 3, it is clear that if 

<3 " w= ^) 8 ° d 

then the overall system reduces to an identity system and hence the reconstruction is exact. 
A few points are noteworthy in this regard. First, Po(z) and P\{z) are assumed to be causal 
stable recursive filters with zeros strictly outside the unit circle and have allpass or near allpass 
characteristics. Second, by definition, the synthesis polyphase filters Qo{z) and Q\(z) will have 
their poles outside the unit circle and their zeros inside the circle. Consequently, they are 
anti-causal and stable. Third, because the image boundaries are finite in length, the circular 
convolution method discussed in [4] is used to handle the analysis and reconstruction at the 
image boundaries in order to avoid an increase in the subband data. It should be noted that 
the recursive polyphase filters are well-behaved and do not attempt to invert stopbands or 
magnify selected regions of the spectrum that might lead to adverse emphasis of coding noise 
[4], Comparisons were given to systems based on linear-phase FIR filters and to those based on 
linear-phase and non-linear-phase IIR filters. The interested reader is directed to this reference 
[4] for further details. 

3 The New Time- Varying IIR Filter Banks 

In this section, we derive a new variant of the IIR filter bank in which the coefficients of the con- 
stituent filters are allowed to vary with time. The coefficients of the polyphase filters Po(z) and 
Pi(z) are selectively changed at some point in the filtering process to another set of polyphase 
filters. This process of changing the filters allows us to improve the coding performance. The 
problem we address in this section is how to reconstruct without distortion once we have switched 
the filter coefficients. By examining the simplified polyphase structure in Figure 3, it is apparent 
that it is sufficient to determine the conditions under which r,[n] = t,-[n] for i = 1, 2. Clearly if we 
can guarantee that r<[n] = <i[n] at all times before, during, and after switching the coefficients, 
then exact reconstruction is also guaranteed. 

To illustrate the analysis and reconstruction problem, consider the graphical illustration 
shown in Figure 4. Assume that the impulse response pi,[n] shown in the figure corresponds to 











that of a causal, stable analysis polyphase filter 


P L {z) = 
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1 + P^Z- 1 + 0%Z- 2 + • • -+ Pk'Z-Ho ' 

At time n = 0, the filter is changed to p/i[n] corresponding to a new polyphase filter 

off + Off* -1 + ^~ Ml 

R ^ Z) 1 + z-1 + P?z~ 2 + ■■■ + /?**“"* • 

The labels L and R signify that the signal or filter coefficient is associated with the left and 
right halves of the time index, respectively. 

3.1 Time Domain Analysis 

The direct form difference equations that implement the polyphase filters are: 
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- m ] - 5j#vz,[n - e\ — oo < n < — 1 
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and 
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Vfl[n] = - m l - 5j9?vjt[n 0 < n< oo (4) 

m=0 /=1 

Here, vl [n] and are the left-sided and right-sided components of the output v[n], i.e. u[n] = 
t)j,[n] + u«[n]. In addition, equation (4) requires the initial conditions uj*[— Ni], v r [-Ni — 1], 
...,«*[—!] in order to evaluate Vfl[0] since 


Mi Ni 

m = 0 1 = 1 


( 5 ) 


Rather than making these initial conditions zero, which effectively discards the recent history 
of the input, we assume the more desirable initial conditions 

Urt[n] = Vi[n] in the range — Ni < n < —1. 

Generalizing the reconstruction procedure discussed in [4] for conventional recursive filter 
banks, we obtain the reconstruction shown in Figure 5 which is based on anti-causal filtering. 
Here #n[n] and gt[n] are impulse responses corresponding to the stable, anti-causal polyphase 
synthesis filters. The corresponding difference equations are 

Mi . Ni 


1 C \ 1 * 

z fl [n-A3=— 1 - + -jp- 'Z2P*i-™ VR [ n + m - ( 6 ) 

a Mi t=1 a Mi m=0 

for 0 < n < oo and 
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x L [n — C] = ~~l~ + + — - m VL[n + m- D] (7) 
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Figure 4: Graphical Illustration of Time- Varying HR Analysis Polyphase Filtering 



Figure 5: Graphical Illustration of Time- Varying HR Synthesis Polyphase Filtering 


for — oo < n < —1, where A, B, C, and D are arbitrary integer shift constants. The reconstruc- 
tion sections, *it[n] and j?r,[n] are right-sided and left-sided components of the reconstructed 
signal z[n], as shown in Figure 5 and have the property that 

j:[n] = x fl [n] + ar L [n]. (8) 

Although the form of the reconstruction is given by equations (3- 7), the unknown parameters 
A, B, C, and D must be determined. 

3.2 z- Transform Domain Analysis 

In this section, we determine the exact reconstruction conditions under this new time-varying 
filter bank paradigm (i.e. find A, B, C, and D ) and show explicitly that reconstruction is exact. 
The time-varying filter bank problem can be analyzed in the z-transform domain in terms of a 
causal and anti-causal unilateral z-transform. We define the causal unilateral z-transform as 

X R (z) = £>[n] (9) 

n=0 

and the anti-causal unilateral z-transform as 

-i 

Ax,(z) = *["]*“"• (10) 

n— — oc 

The sum of Xr(z) and Xl{z) is identical to the bilateral z-transform of x[n\. These causal and 
anti-causal transforms have shift properties associated with them that may be described in the 
following way. Suppose that no is an integer and u[n] is the unit step function. Then, for a time 






shift of no in z[n] 

to[n] = x[n - no] 

the causal unilateral z-transform W R (z) is 


n 0 — 1 


— n 0 — 1 


W R (z) = z~ n °X R (z) + *l £ ~ noKMno - 1] - 2 '"° 22 *M*"M-no - !]• 

i=o i=o 

Similarly, for the anti-causal z-transform, the same time shift 

tu[n] = x[n - no] 


( 11 ) 


results in the anti-causal transform 

no 


— fio — 1 


Wl{z) = ^ n °X L (z) - z n °22 x [— £ \ ziu [ n o - 1] + z~ n ° 22 *M 2- * u [ _n o - 1]. (12) 

/=i /=o 

This shift property will be used to take the transforms of the analysis equations (3) and (4) and 
the synthesis equations (6) and (7). A derivation of this shift property is given in appendix A. 

Using the causal and anti-causal z-transforms and their corresponding shift properties, we 
can take the z-transforms of the analysis and synthesis equations. The transform properties 
allow us to model the initial conditions explicitly for the point at which we switch the filter 
coefficients. In this analysis, we assume that we switch the filter coefficients at n = 0. 

For convenience in taking the transform, we express equation (4) as 


Mi 

22 /3 * VR [ n - *1 = £ a m x l n -™] o < n < oo 

1=0 m = 0 


(13) 


(where /3q = 1) and compute its z-transform using the definition (9) and shift property (11). 
We obtain 

Ni Ni l-l Mi Mi m— 1 

22^z- t V R (z) + 22^T / v ^ k ~^ z ~ k =12 Q ^XB( z )z-’ n +22^ J Z x 4 £ -rn)z- t . (14) 

1=0 1=1 k =0 m = 0 m = 1 /=0 

Next we can write the right-sided synthesis equation (equation (6)) as 
Mi Ni 

22 XR ^ n+t ~ A ^ a Mi-i = J2^i- m M n + m ~ B ] 

1=0 m =0 


and take its z-transform to obtain 


A-t-l 


Mi Mi 

22°‘Mi-i x r{ z ) z1 ~ A + 22 a Mi-i( £ x R[™-A + e\z~ m u[A-e- 1] 

r=0 1=0 m = 0 


— z 


-yi+r-i 

- A+< 22 ** m « m +*- id 

m=0 


Ni 

I 

m =0 


Ni 


B—m—1 


= £ 0*i-Mz)z m - B +22/*Ni-m( 22 + 1 ] 


m =0 


/=0 


(15) 


(16) 


—B+m—l 

£ t>*[<M-B + m-l]). 

1=0 

This is the z-transform of the synthesis equation for the right-sided sequence Xfl[n], If the 
right and left sides of this equation can be shown to be identical, then we have shown exact 
reconstruction for zjj[n], To explicitly show that the equality holds, we must express Vr(z) in 
terms of Xr(z), the initial conditions, and the filter coefficients. Using equation (14), we can 
remove the Vr(z) term by substituting 

Ml Ml 771 — 1 -Wj /—I 

- i\z~ k ) (17) 

m= 0 m=l 1=0 1=1 k = 0 


for 


N, 




m—B 


m=0 


in equation (16). After the substitution, equation (16) becomes 
Mi Mi A-t-1 

Y a Mi-t x R( z ) zt ~ A + Y a Mi-i Y XR[*n~A + e\z- m u[A-l- 1 ] 


1=0 


1=0 


m=0 


( 1 ) ( 2 ) 

Mi -A+l-1 

-Y a Mi-t*~ A+t Y **M«[-^+*-i] 

1-0 m= 0 

' — ' 

( 3 ) 


Mi M\ tii—1 Afj /—I 

= Z Nl ~ B Y a m^fl( Z ) 2 " m + Z Nl ~ B Y°*Y XR [ £ ~ m ^ Z ~ l ~ zNl ~ B YP*Y VR { k ~ 4 z ~ k 

m=0 m=l /=0 t= 1 k = 0 


( 4 ) 

B — m—1 


( 5 ) 


( 6 ) 


N i B—m—1 Ni —B+m—l 

+ YP*i~™ Y '>R{t-B + m}z- l u[B-m-\]-Yl }B i-mZ B - m Y v R [£\u[-B + m - 1] 

e=0 ^ m=0 t= 0 


m = 0 


( 7 ) 


( 8 ) 


(18) 

By examining this equation, it becomes apparent that we can group like terms and that there 
must be equivalence within these groups in order for the equality to hold. Consider, first, terms 
6, 7, and 8, which are the only terms containing components. The sum of these terms 

must be zero, i.e. 

Ni l- 1 Ni B-m-1 

-z N '- B Y,0*Y VR l k - + Y^i-m Y VR[l-B + m)z- l u{B-m- 1] 

^ <_i k= o ^ ^=0 r=o 


(6) 




Ni —B+m—l 

-Y$i-™ zB ~ m Y v R m-B + m-l] = 0 (19) 

m=0 <=0 

v 7 

(8) 


At initial glance, it is apparent that cancellation can only occur between either terms 6 and 7 
or 7 and 8 because these are the only combinations in which the signs are opposite. A closer 
examination shows that terms 7 and 8 are disjoint in the variable m due to the u[B — m — 1] 
and u[—B + m — 1] terms. Thus, it is sufficient to focus on terms 6 and 7. In term 6, if we let 
£ = N\ — m we obtain 

JV,-1 

m=0 i=0 

which cancels with term 7 when B = N\. Moreover, when B = N\, term 8 disappears because 
t »[— B 4- m — 1] goes to zero. 

We can now group terms 1 and 4 together in equation (18) because only they contain the 
Xr{z) components. Letting £ = Mi — m in term 1, we obtain 

Y<**X R (z)z M '- m ~ A . 

mz rO 

This becomes identical to term 4 when A = Mj, resulting in cancellation. 

The only remaining terms are 2, 3, and 5. Immediately we observe that term 3 is zero. Note 
that because A = Afj, the u [— A + £ — 1] forces term 3 to be zero. Thus all that is left is to show 
that terms 2 and 5 cancel. Letting £ — — n in term 2, we obtain 

Mi n — 1 Mi n — 1 

n]z~ m «[n - 1] = Y a * n l z ~ m 

n=0 m=0 n=l m=0 

which is equivalent to term 5. Thus we see that when A = —Mi and B — —Ni , we can exactly 
reconstruct 

What remains is to prove that we can reconstruct Xi[n], The same procedure can be applied 
to the left-sided equations. In particular, the analysis equation (3) can be written as 

No Mo 

XX V£ [ n = a ^ n ~ ( 2 °) 

1=0 m= 0 

where = 1. Its corresponding z-transform is 

No No e Mo Mo ro 

z ~ l yL(z)~Y0i z ~ t 'Y VL ^~ k ^ zl ‘ u ^~ x i= Ys a ™ z - mXL ( z }~Yj a ™ z ~ m '52 XL {~ e \ zl ' i \ rn - 1 } 

1=0 1=0 i = l m = 0 m= 0 t= 1 

(21) 

The z-transform for the left-sided synthesis equation (7) is 

Mo Mo C-t t-C-1 

-tXL( z ) zl ~ C - ^2 a Mo-l zl ~ C (Y^ XL ^~ k ^ zku ^ C ~ 1 ~ ^ + YL x L[k]z- k u[£ - C - 1]) = 

/=0 t= 0 k=l k= 0 

No N 0 D-m —D+l—1 

- m zm - D V L {z) - Yrto-m zm - D ( Y ~ ™ - 1} A Y ~D- 1 ]) ( 22 ) 

m=0 m=0 l—\ /=0 


Applying equation (21) to equation (22), we can remove the Vi(z) components and we obtain 


M 0 Mo C-t 

Y^ a M 0 -tXUz) zl ' C - ~ 1 - !] 

i = o r=o t=i 

l-C - 1 

+ £ x L [k)z- k u[t-C-l]) = 

k = 0 

No l M 0 

z N '- D (^z- t '£vL[-t]z k u[i- 1 ]+ J2^~ m X L (z) 


t = 0 i=l 
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( T: -m - 1]+ y v L [C]z~ l u[m- D - 1]) (23) 

/=i r=o 

Again we can equate like terms and solve for the unknowns. Doing so, we find C = Mo and 
D = No, and consequently prove that £/,[«] may be exactly reconstructed. 


4 Results 


The time- varying filter bank reconstruction problem shares similarities with conventional inverse 
polyphase reconstruction. However, the interactions at the filter coefficient switching points 
make the analysis very different for the time-varying case. Analysis and interpretation of time- 
varying filter banks is complicated by the fact that the tools we would like to use, such as the 
DTFT and 2 -transform, are inherently time-invariant. By couching the problem in terms of left 
and right signal components and the initial conditions associated with them, we have proved 
that maximally decimated time- varying filter bank analysis/synthesis systems can be exactly 
reconstructing. We determined that the unknown parameters in the analysis equations (3) and 
(4) and synthesis equations (6) and (7) are 

A = Mi 
B = Ni 
C = M 0 
D = N 0 

(24) 


and that the resulting reconstruction synthesis equations are 

M t , Ni 


1 ' % 1 ^ ^ 
x R [n - M x ] = g- + + X IpNi-mMn + m-Ni] 
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(25) 


‘Mi t=1 
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for 0 < n < oo with t>jt[n] = vi[ri\ in the range —N\ < n < — 1 and 

^ Afo ^ No 

x L [n — Mo] = — J2 a Mo-t x ^ n + e -Mo] + -ZIT Yl0No-mVL[n + m- N 0 ] (26) 

~ a M 0 r = i “Mo m=0 

for — oo < n < —1. 

4.1 Constraints on Switching in the System 

This 2 -domain analysis and the resulting synthesis equations reveal some interesting and impor- 
tant properties of the analysis/synthesis system. To begin, observe that the derivation leading 
to equation (25) proves that the samples of the input x[n] in the range —Mi < n < oo can 
be reconstructed using equation (25), where z[n] = z«[n]. Similarly the samples of x[n] in the 
range — oo < n < —Mo can be reconstructed using equation (26), where x[n] = n[n] in this 
region. 

In order to reconstruct all sample values using equations (25) and (26), constraints upon the 
variables Mo and Mi must be imposed. This observation can be drawn from Figure (6) with 
the following constraint being placed upon Mo and Mi : 

Mi + 1 > Mo- 

If this constraint is not met, the values in between Mo and M\ can not be recovered from 
the reconstruction equations (25) and (26) as shown in Figure (6). There exists a gap of 
unrecoverable samples in between Mo and M\ . 

Equations (25) and (26) also show that by supplying the original sample values in the 
transition regions explicitly, the constraint on the separation between switching points can be 
removed completely. As we saw, switching coefficients at an interval less than or equal to Mi 
precludes the recovery of all samples in the region —Mi < n < 0. If these missing samples were 
supplied externally, then reconstruction can proceed without difficulty. Unfortunately, such an 
approach results in a data increase, i.e. the analysis/synthesis system is no longer maximally 
decimated. 

For convenience of being able to use the conventional unilateral 2 -transform, the analysis 
has been given under the assumption that the switching of filters occurs at n = 0. But nothing 
is changed by switching the coefficients at an arbitrary time, n = no, except for a time shift. 
All the equations derived are applicable to a system with arbitrary time shift at n = n 0 . Based 
on the analysis, it should be clear that these results hold independent of when the filters are 
switched. 

However, the equations reveal limits on how frequent these switching of the filter coefficients 
can occur within a given interval and on the numerator/denominator orders of the polyphase 
filters. To illustrate this, consider the example in Figure (6). Equation (25) reconstructs valid 
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Figure 6: Reconstructing Samples using Equations 25 and 26 Illustrating Constraint on Mo and 
Mi 


samples of Zfl[n] = *x,[n] for n = -1,-2, —Mi. These samples are used as initial conditions 
in equation (26). However, if another switch of coefficients occurs at time n = K where —Mi < 
K < 0, only K valid samples can be reconstructed from equation (25), the remaining samples 
in the interval will be erroneous due to the switch in coefficients. Thus, the analysis proves 
that exact reconstruction can be achieved when the filters are changed at arbitrary intervals 
provided that the switch points are separated by a number of samples equal to or greater than 
the numerator order of the associated analysis polyphase filter. Assuming and nj+i are 
consecutive switch point indices in n, the criterion for exact reconstruction is 

n,+i - n, > Mi. 

We hasten to point out that this is an extremely mild restriction in practical applications. This 
is because the orders of practical recursive polyphase filters tend to be very low. Recursive filters 
which are first-order allpass polyphase filters can have magnitude characteristics comparable to 
a 24- or 32-tap QMF. Such filters have extremely good magnitude characteristics and generally 
provide more than enough stopband rejection with sufficiently narrow transition bands. Because 
the polyphase filter order is just one in this case, the polyphase filters can be switched every 
other sample if needed. 

4.2 Stability Issue 

Another issue to address is stability. It is well known that if the filter coefficients are allowed 
to change arbitrarily, the system is not guaranteed to be stable. In our case, however, this is 
not a problem. In the formulation for the recursive time-varying analysis-synthesis system, the 
general problem is decomposed into piecewise constant coefficient system components, each of 
which is stable^ince the switch points are separated by Mi samples, each region (after a switch) 
that is reconstructed by equation (26) can be represented uniquely by a unilateral z-transform 
term. In fact, all reconstructed regions that lie between switch points can be represented by 
separate unilateral z-transforms. 




To illustrate this, consider the first-order allpass polyphase filter all-pass polyphase filters 

_ . . a — z -1 , . 

P o( z ) = Y - az -~ ' N <1 * 

which we assume to be a filter employed over a finite interval. The corresponding difference 
equation is 

y 0 [n] = ay 0 [n - 1] + ax 0 [n) - z 0 [n - 1] 


and is used to generate all the samples in the interval. The initial conditions for the equation 
are derived from the output in the previous interval and in general will not be zero. If we take 
unilateral z-transforms on the above equation, we obtain the transform-domain equation 


Yo( z ) = [*«(*) + 


at/o[-l] - *o[-l]' 


- 2 -i 


where t/ 0 [ — 1] is the initial condition. Since, the value of a is less than one, the equation output 
will always be bounded. 

The stability is thus completely determined by the locations of the poles of the corresponding 
transform. That is to say, for stability the poles associated with the analysis polyphase filters 
should be inside the unit circle and the zeros of the analysis polyphase filters (which become 
the poles of the synthesis polyphase filters) should be outside the circle. 


5 Experiments and Performance Analysis 


The analysis presented here used the simple case of a two-band polyphase system to develop the 
reconstruction equation derivation. Since the result applies to the actual polyphase and inverse 
polyphase filtering operations, banks with complex- valued channels [11] and an infinite variety 
of tree structures composed of two-band and/or parallel form filter banks. 

Finally, as alluded to in the introduction, time- varying filter banks can improve the subjective 
and objective performance because of the input dependency of the error characteristics. This is 
clearly the motivation for the development of time- varying filter banks. It is also a measure by 
which one can judge its utility. 

As an initial gauge of performance, we compare a conventional two-band subband coder 
and a two-band subband coder with time-varying filter coefficients. A finite length sequence 
was used as the input — shown in Figure 7 as the solid line. This input is a section of a row 
taken from the 256 x 256 test image “Lena.” In each case, the lowpass channel was coded 
with 4-bit uniform quantization while the highpass channel was coded with 2 bits. Notice that 
the conventional subband coder displays ringing (or overshoot-undershoot) distortion at the 
discontinuities as indicated by the dashed line in Figure 7. However, the other regions of the 
row are represented well. The magnitude and step response characteristics for Ho(z) used in 
the conventional subband coder are shown in Figure 8. 



Figure 8: Magnitude and step response for the lowpass filter II 0 (z). 


The new subband coder with time-varying coefficients employs two filter sets, one with good 
magnitude response characteristics (i.e. Ho(z) shown in Figure 8) and the other with good 
step response characteristics. The magnitude and step response characteristics for the latter 
filter are shown in Figure 9. In the regions of large discontinuity, the filter set with good step 
response characteristics is used. In other regions, the filter set with good magnitude response 
characteristics is used. Figure 10 shows the improvement that results for this example. 

6 Remarks 

Such an approach involves having to convey switching point information to the receiver for 
reconstruction. This can be done simply by sending side information as was the case in the 
example or potentially through a procedure based on feedback, the latter of which is presently 
being studied. For natural images, the amount of side information is negligible. Thus, even 
with the transmission of side information, this new approach seems attractive. 
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Appendix: Derivation of Shift Property 


The unilateral z-transform is defined as 

X R (z) = f 

n=0 

If z[n] is shifted by a positive integer time shift TV, and the transform computed, we obtain 

f>[n-TV]z'". 

n=0 

Substituting n = ifc + TV, results in 

CO CO —1 

53 *[*]*"*-" = z-*]jr *[*>-* + z-* *m*~‘ 

k=-N k=0 k=-N 

= X r (z)z~ n + z~ N £ *[fc]z- 1 . 

k=-N 

If we let k = TV — i and substitute, we obtain 

N - 1 

X*(z)z-* + ]z' £ . 

1=0 

Because this is valid only for integer N > 0 we can append the term u[7V — I], For negative 
shifts, TV < 0, the transform becomes 

co oo co — TV— 1 

£*[«-7V]z- = £ r[i]r‘-« = r£#-‘- 2 -^ £ 

n=0 k=-N k=0 k= 0 

-N- 1 

= X fl (z)z- N - Z-* ^ x[fc]z~*u[— TV — 1] 

k= 0 

where u[— TV — 1] restricts TV to be less than zero. Combining these results yields 

N-l -N - 1 

x[n - TV] <=> z~ N X R {z) + z~ N ~ N]z-*u[N - 1] - z~ N ]T z(fc]z-*u[-TV - 1]. 

1=0 k=0 

For the anti-causal case, the derivation is similar. The transform of the time shifted input is 

-l 

^ z[n-TV]z _n . 

n=— OO 

Assuming TV > 0 and substituting n — k + TV, results in 

£ zfjfcjz - * - * = z~ N £ x[k]z- k - Z ~ N £ x[k]z~ k 

ifc=— oo kzz—oo k=—N 

N 

= X l (z)z~ n - z~ N ^2x[-k ]z k 
k = 1 


For N < 0, the transform becomes 


E 

n =—oo 


-N - 1 


-1 


-JV-l 


53 *w *"*' JV = *~ N E *[*k*+* _n E 


k = — oo 


k=0 


Combining both results we obtain 

N -N-l 

x[n -N]<=* z~ n X l {z ) - - 1] + z~ N £ *[t\zr l v[-N - 1]. (27) 

r=i r=o 


